Hidden spin liquid in an antiferromagnet: Applications to FeCrAs 
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The recently studied material FeCrAs exhibits a surprising combination of experimental sig- 
natures, with metallic, Fermi liquid like specific heat but resistivity showing strong non-metallic 
character. The Cr sublattice posseses local magnetic moments, in the form of stacked (distorted) 
Kagome lattices. Despite the high degree of magnetic frustration, anti- ferromagnetic order develops 
below Tjv ~ VSiK suggesting the non-magnetic Fe sublattice may play a role in stabilizing the 
ordering. From the material properties we propose a microscopic Hamiltonian for the low energy 
degrees of freedom, including the non-magnetic Fe sublattice, and study its properties using slave- 
rotor mean field theory. Using this approach we find a spin liquid phase on the Fe sublattice, which 
survives even in the presence of the magnetic Cr sublattice. Finally, we suggest that the features 
of FeCrAs can be qualitatively explained by critical fluctuations in the non-magnetic sublattice Fe 
due to proximity to a metal-insulator transition. 



I. INTRODUCTION 



The ubiquity of Landau's Fcrnii liquid is a testa- 
ment to universality in the solid state. As such, de- 
partures from these classic experimental signatures in 
metallic systems act as a guide to novel and interesting 
physics. Non-Fermi liquid behaviour appears in many 
strongly correlated materials such as un con ventional 
superconductorgilEl, heavy fermion materials^ and near 
quantum phase transitions -^. Some routes to realize this 
behaviour include coupling itinerant electronic systems 
to localized magnetic moments and through intermedi- 
ate to strong electron-electron interactions. These mech- 
anisms can give rise to characteristics and experimen- 
tal signatures that do not fit neatly in the Fermi liquid 
paradigm. 

A recently re-examined compound, FeCrAgS'ISl^ pro- 
vides a direct example of a material that does not fit 
completely within Fermi liquid theory and combines as- 
pects of the mechanisms discussed above. The unit cell of 
FeCrAs shown in Fig. [l] shows Cr and Fe form alternat- 
ing two dimensional lattices along what we will denote 
the c axis. Cr forms layers with the structure of a dis- 
torted Kagome lattice where the Cr — Cr distances are 
approximately constant. The Fe layers have a more com- 
plicated structure, forming a triangular lattice of three 
atom units which we will call trimers as shown in Fig. [2] 
The As is interspersed throughout both the Fe and Cr 
layers, as well as in between. Most of the known experi- 
mental data is nicely presented in Wu et aP^, which we 
summarize below. 

The specific heat exhibits Fermi liquid behaviour at 
low temperatures, i.e. C ~ 7T where the slope 7 is sam- 
ple dependenlflH. The measured linear range is roughly 
T ~ iK — VdK. Resistivity measurements show insulat- 
ing behaviour at low and high temperatures. In plane 
(a6) and out of plane (c) resistivities are of the same or- 
der over the entire temperature range considered. The 
resistivity monotonically decreases (that is dp/dT < 0) 
as T is raised from 5K up to ~ 800K except for a small 
peak in the c axis resistivity around T ^ 125K. A low 








(a) View along c axis 



(b) Tilted view 



FIG. 1: (Color Online) FeCrAs unit cell viewed along the 
c axis (left) and tilted away by 70° (right). Fe atoms are 
indicated in green, Cr in blue and As in yellow. 





(a) Trimer Lattice 



(b) Relation to Cr sublattice 



FIG. 2: (Color Onhne) (a) The Fe sublattice with trimers 
shown explicitly, (b) The local environment of the trimer, 
with respect to the Cr sublattice. 



temperature power law p ^ po — AT" is observed for 
T ~ 80mK — 5K in both ab planes and c axis resistivity 
with a - 0.6 - 0.7 

There is a peak in the susceptibility at T/v ~ 125K 
indicating a magnetic transition with a lack of hystere- 
sis pointing to antiferromagnetic ordering. Below T^v 
the susceptibility is anisotropic, differing between the ab 
plane and the c axis. Elastic neutron scattering^^ done 
deep in the magnetic phase, at T = 2.8K, is consistent 
with the anti-ferromagnetic order inferred from the sus- 
ceptibility, signaling an ordering vector at Q — ( | , | , 0) , 
indicating ferromagnetic (stacked) order along the c-axis, 
but with a tripled unit cell in the ab plane. While 
measurements of the specific heat give a result consis- 



tent with a metallic Fermi liquid, transport is unusual 
and deviates strongly from the classic Fermi liquid result 
for metals, while being distinct from the expected result 
for strong insulators. Furthermore, this material has a 
frustrated magnetic sublattice that nonetheless orders at 
low temperatures, while the remain ing su blattices show 
either small or no magnetic moment ^^ I The magnetic 
sublattice takes the form of a distorted Kagome lattice, 
where even classical Heisenb erg m odels fail to order mag- 
netically for both the stackecP^lMl and purely two dimen- 
sional casesi^H^. Very few experiments have been carried 
out on FeCrAs, so theoretical models are not completely 
restricted. Regardless, there are a number of questions 
that need to be addressed, such as the nature of the stabi- 
lization of magnetic order, the cause of the very different 
thermodynamic and transport signals and the role of the 
non-magnetic sublattice. 



The nature of the magnetic order has been recently 
addressed by Redpath et aP^l^ where a minimal model 
was proposed which suffices to explain the experimentally 
observed stabilization of a particular magnetic ordering 
vector. However, transport and thermodynamic behav- 
ior remains to be explained along with the role of the 
non-magnetic, Fe sublattice. In this paper, we elaborate 
a microscopic route to an effective model for the com- 
pound FeCrAs, taking into account the Fe sublattice, and 
present a scenario to address the incongruities between 
the conflicting metallic, Fermi liquid specific heat and 
insulating like transport signals. This model consists of 
interacting electrons of the non-magnetic sublattice cou- 
pled to magnetic moments of the magnetic, Cr sublattice. 
Here we do not address the detailed nature of the mo- 
ments themselves, treating them classically^*', with the 
non-Fermi liquid physics arising from strong charge fiuc- 
tuations occuring at intermediate Hubbard coupling in 
the Fe sublattice. In this picture we are excluding any 
Kondo physics, a view supported by the experimental re- 
sults. Our emphasis is on the interplay between strong 
charge fluctuations near the metal-insulator transition on 
the Fe sublattice and the magnetic order of the Cr sub- 
lattice. For this we turn to the slave-rotor method which 
allows access to the intermediate coupling regime and 
metal-insulator transition. 



The structure of the paper is as follows: in Section |TT] 
we present an argument to pass from the atomic limit 
through to an effective model of the electronic degrees of 
freedom in FeCrAs. In Section Hi Dl we discuss the local- 
ized moments and magnetic interactions and we present 
the effective Hamiltonian relevant for FeCrAs. We pro- 
ceed to review the slave-rotor method in Section IIIII and 
the assumptions and implementation of our mean field 
theory in Section |IV| In Section |V] we comment on the 
application of our results to FeCrAs and summarize con- 
clusions in Section IVTl 





(a) Tetrahedral (b) Pyramidal (c) Trimer 

FIG. 3: (Color Online) The local environments of the Fe (a), 
Cr (b), and trimer (c). The Fe atom is tetrahedrally coordi- 
nated and the Cr atom is approximately octahedrally coor- 
dinated. The trimer is surrounded by three tetrahedra that 
share a common axis. 
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(a) Effect of the 
tetrahedral crystal field 



if.S^^^ 



(b) Low spin state for 
Fe3+ 



FIG. 4: The effect of the crystal field on the Fe^+ ion. The 
local symmetry about this site is tetrahedral, that is the group 



II. EFFECTIVE HAMILTONIAN 

A. Local environments and spin states 

Considering the common oxidation states of Fe we will 
take (in the atomic limit) Fe'^^ as a starting point. This 
leaves the valence configuration being 3d for Fe^^. Fol- 
lowing the experimental data, we assume that the Fe 
atoms are in a low spin state, due to the lack of a de- 
tectable magnetic momenlfi^l. The tetrahedral arrange- 
ment of the As atoms about each Fe atom, shown in Fig. 
[3] (a) , implies a crystal field producing a splitting of the 
Fe d levels, shown in Fig. [4] (a), with a pair of low- lying 
e levels and a three-fold degenerate set of ^2 levels, sep- 
arated by the crystal field gap. Assuming the low spin 
case is relevant the two e states are filled and there is a 
single electron in the ^2 triplet shown in Fig. H (b) . 

The case of Cr is more complicated, as the experiments 
do not single out a more probable spin state. For Cr we 
will take a ionic charge of Cr , giving a valence con- 
figuration of Sd"*. A distorted octahedral environment, 
shown in Fig. ^ (b) reduces the symmetry about this 
site to D4, with crystal field splittings shown in Fig. [5] 
(a). The low spin state(shown in Fig. [s] (b)) can still 
yield a S" = 1 spin moment, while the high spin state has 
an 5 = 2 moment. In the low spin case hopping onto the 
Cr^+ will be suppressed by the orbital repulsion, while 



ai 



'A 



D4. 



I = 2 

t2g 



61 



ai 
fei 



A' 



\^n. 



(a)Effect of the tetragonal crystal 
field 



(b)Low spin state for 
Cr2+ 



FIG. 5: The effect of crystal fields on the Cr^+ ion. The Cr 
ion is biased towards one direction of the octahedron. This 
reduces the symmetry to the tetragonal group D4. 



for the high spin case the crystal field energy will give a 
further supression. 

We would like to emphasize that the arguments and 
models presented here are under-constrained by both ex- 
perimental results and first principles electronic struc- 
ture calculations^]. Due to this limitation the precise 
details could fail quantitatively, but the subsequent ef- 
fective Hamiltonian appears to be robust to a variety of 
ionic configuration changes in the underlying model. For 
example, the specific oxidation state we use for the Cr 
will be irrelevant to our final discussion, as only the lo- 
calized character is needed to capture the gross magnetic 
featurea^Sl. With this in mind below we present a possible 
route from the microscopic Hamiltonian to the effective 
model. 



B. Iron-Chromium interactions 

Since the Fe-Cr distance is considerably smaller than 
the direct Fe-Fe distances outside the trimers and the 
indirect Fe-As-Cr distance, we will consider interactions 
induced by Fe-Cr-Fe hopping paths. First we define a 
simple model for an isolated Cr atom, assuming a low 
spin state as shown in Fig. [s] (b) . The local Hamiltonian 
for the Cr e doublet assumes the natural form 



Fe site i. Integrating out all of the high energy degrees of 
freedom on the Cr atom, the ground state for sufficiently 
large J is given by the doubly occupied, S = 1 triplet 
states which we denote by \tja) with a — Q,zt. Noting 
that _ffci._Fc does not connect triplet to triplet states, as 
it changes the electron number on the Cr atom we can 
formulate the effective Hamiltonian following 



Hr 



PtHci-Fc {Eq — Hcr) 



-"Cr- 



Fo-ft: 



(3) 



where Eq is the triplet energy and Pt = X^ia Kja) ('•jal 
projects onto the triplet subspace. Carrying out the ex- 
pansion of the effective Hamiltonian, denoting the spin 
operators on the Fe and Cr atoms as Si and i^ai respec- 
tively, our full effective Hamiltonian is given by 
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(5) 



For i = I this represents an anti-ferromagnetic exchange 
between the Fe and Cr atoms, while i 7^ I presents a 
spin-dependent hopping term, non-zero for both intra- 
trimer and inter-trimer hopping paths. Based on overlap 
of the orbital wavef unctions, we assume that the hopping 
J2a l^fiP is dominant for i = / and independent of the 
orbitals, leading to the effective exchange Hamiltonian 

^eT'^ = E^fEi^Si')4(-5'.)rf. (6) 

ij \ a I 

up to a shift of the chemical potential. For future use, 
we will denote the effective exchange as 



'^^ = ^Ei^si'>o- 



C. Trimer Approximation for Fe 



(7) 



TJcr = A E '^j" + ^ E "j"t"ja; - -^ E ( E '^j" ) ' 

ja ja j \ a / 

(1) 

where Ujacr — ^jaa^jacr IS thc number operator for the 
state in the doublet a. We have denoted the atomic 
potential as A, the intra-orbital repulsion as U and the 
Hund's coupling as J. Allowing for hopping between Fe 
and Cr, we include the term 



ii^Cr-Fo = - E *S4<Teiao- + h-C, 



(2) 



where i runs over the orbitals on sites connected to the Cr 
atom at site j and d is the annihilation operator on the 



From the inter-atomic distances we expect that the Fe- 
Fe atoms in the trimer structure are tightly coupled. We 
can take advantage of this by grouping the degrees of free- 
dom in the trimer into a single unit and formulating the 
rest of our theory in terms of these variables. This pro- 
cess is similar to the treatment of the pairs of molecules 
as the effective degrees of freedom, frequently employed 
in studies of organic superconductors'^ . This assumes the 
energy scales for interactions and inter-trimer hopping 
matrix elements are small relative to the intra-trimer 
hoppings. In addition, we need that the tetrahedral crys- 
tal field splitting be much larger than our intra-orbital 
interactions (this is the low spin assumption) and both 
the intra- and inter-trimer hopping elements. 



: t^^, 



of the three sites is reahzed. Ignoring ah of the other 
matrix elements in the complete hopping matrix, we find 
(see Appendix A) 
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(a) Spectrum when 

ixz,yz ^ ixy ,xy 



(b) Spectrum when 

ixy,xy ^ T^xz^yz 



FIG. 6: The two cases that give rise to a half-filled band 

with t 

character or with (b) d+ character. 



xy.xy ~ txz,yz- We Can either have a band with (a) dxy 



With this in mind, let us find the low energy degrees of 
freedom of a trimer, keeping only the three site cluster. 
Considering the C3 symmetry, there are four possiblities 
for the degeneracies of the two lowest levels, leading to 
either a half or quarter filled band as the relevant states 
(assuming all gaps are large compared to the relevant 
energy scales) . Motivated by the Slater-Koster argument 
in Appendix A, we start with a model for the xy, xz and 
yz orbitals of the trimer: 
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(8) 



txz,vz ^ ^ {7ts - 16<^ + 9tcr) 



''xy.xy ^ „, (.4915 



I2t^ + 3t„) 



(11) 
(12) 



A more complete model would include mixing between 
all three orbitals but the qualitative picture should re- 
main the same. The naive choice of t^ > i^ > io- seems 
to select the case of t^y^xy > txz,yz and thus the single 
occupied band has d+ character (shown in Fig. [6]). Note 
the gap between the occupied level and the half-filled 
level is given by 



At 
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rnncr '^\^xy,xy 



^xz^yz 
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-4i, 



3t.|, (13) 



and must be a fairly large to ensure that the trimer ap- 
proximation is valid. We emphasize that regardless of the 
exact details of the these hoppings the qualitative picture 
will remain, as long as this half- filling is found such as for 
Fe'^+ . Other oxidation states will give a different relevant 
molecular orbital, possibly with extra orbital degeneracy, 
and could change the picture presented here. 



where d] ^ { d. 
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-i.xz -^^yz -i^xyj and i, j = 1,2,3 are the 
sites in the trimer. The xy part is just a simple three 
site chain with a ground state at —2txy^xy and pair of 
excited levels at txii.xv For the xz and yz orbitals we 



xy,xy- 



can diagonalize -fftrimer by changing the basis: 



1 



di,+ — —7= [dxZ + dyz) 

v2 

^i.~ T= \(^xz ^yz) 

v2 



(9) 
(10) 



This gives a pair of decoupled three site chains with hop- 
Hxz,yz and thus energy levels of ±2txz-yz and 
< txz,yz < 2txy,xy a single half- 



'^T'xz^yz- VV nen 2^xy,xy 



filled level in the trimer, symmetric under permutations 



D. Effective Hamiltonian and magnetic interactions 

The gross magnetic structure of FeCrAs is captured by 
a classical model of localized moments interacting with 
the Fe sublattice via a simple exchangers. The utility 
of the classical model for the moments is also supported 
by the fact that Kondo-like signatures are absent in the 
experimental data. Futhermore, due to the uniform char- 
acter of the relevant trimer states and the assumption of 
equal exchanges between Cr and different Fe orbitals, the 
effective exchange between the Cr and the trimers will 
be equal and between all of the Cr in the surrounding 
hexagons in the layers above and below. We thus take 
our Hamiltonian to have the form 



Jk 



H = -t J2 E 4.^^- - *' E E ^l^^- + f^ E "»t"4 + ^ E (^I'^^O ■Sa + jHY.Sa- s,, (m) 



{ij)Gab cr 



{ij)ec <y 



i.aeOi 
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where the low energy degrees of freedom on the trimers afforded hoppings t and t' which originate fronr direct 
are denoted using the operators c|^ , Cta and the classical overlaps between trimers (or indirectly via Cr or As) in 
Cr spins are denoted as Sa- The trimer sublattice is ^^^ «^ Planes and along the c axis respectively. Here we 



denote sites on the Fe sublattice as i and j, while sites on 
the Cr sublattice are denoted as a and b. The notation 
a E Oi indicates that the sum on the Cr sublattice is over 
sites in the hexagon that surround the trimer at i on the 
Fe sublattice. Furthermore (ij) € ab, c denotes bonds 
in the ab or c planes respectively. The interactions are 
given as follows: U is the intra-trimer repulsion, inherited 
from the Fe atoms, Jk is the Fe — Cr exchange and Jh 
is the Cr — Cr exchange. In this Hamiltonian the Cr 
spins appear as an effective magnetic field for the trimers, 
which we will denote as h^ = X)agO ^o-- 

In the large U hmit, where f/ > t and t'^/U, {t'Y /U < 
Jk,Jh, this model should reproduce the model studied 
in Ref |20j . and thus their classical results provide a use- 
ful point of comparison for our large U /t behaviour. To 
attack the intermediate U /t regime we will use a slave- 
rotor approach, reviewed in the following section. 



III. SLAVE ROTORS 

For completeness and standardization of notation, we 
review the general properties of two dimensional rotors 
and follow with a discussion of the slave particle repre- 
sentation that bears their name^^l. 

A rotor in two dimensions is an object that possesses 
only angular momentum, prototypically of the form H ex 
L^ where L is the angular momentum operator about 
some axis, say the z axis. We classify states by the 
eigenstates of the L operator, L\n) = n \n) where n is 
an integer. Raising and lowering operators are defined as 



L/t 



1) 



U\n) 



1) 



where C/ is a unitary operator. From this definition it 
is simple to show that L and U satisfy the commutation 
relations. 

Since U is unitary it can be written as U — exp(— 10) 
where 9^ — 9 and one can show that this implies the 
canonical commutation relation [9,L] = i, showing that 
L and 9 are canonically conjugate variables. 

To use these rotors as a slave-particle we associate the 
local electron basis with the product of the states of slave 
fermion and the states of an 0{2) rotor, 
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The slave fermion is called a spinon and spinon states 
are denoted by an / subscript. The slave-rotor will be 
referred to as a rotor and a 9 subscript will be used to 
denote rotor states. The natural interpretation is to have 
the spinon to be neutral and the rotor to carry the charge 



of the electron, thus explicitly separating the spin and 
charge degrees of freedom. Having expanded our local 
Hilbert space, a constraint is required to remove the un- 
physical states. The four physical states above are char- 
acterized by 



E^ 



,/ 



1, 



where Li is the rotor angular momentum operator and 
nj^ is the spinon number operator. This is the Hilbert 
space constraint. The electron operators can therefore 
be expressed as 



where exp(— i^^) is the rotor lowering operator and fi^-, 
//g. are the fermionic spinon operators. Using this repre- 
sentation the electronic Hamiltonian for the trimers ( 14 ) 
is written as 
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The Hubbard term is now a kinetic term for the rotors, 
so the complexity of the Hubbard interaction has been 
moved to the hopping term and the constraint. We note 
that the Fe — Cr coupling term only involves the spinon 
degrees of freedom. 



IV. MEAN FIELD THEORY 

We approach this problem using mean field theory. 
For simplicity we take the perspective of Florens and 
Georges,'^ first decoupling the hopping term into 
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ij, 



where we have introduced the mean fields Xij = 
5 EaiflM and By = (e-*(«'-«^)). Note that we have 
assumed that Xij is independent of spin. To handle the 
constraint we treat it both on average in space and on 
average in our states. For the case of half-filling this leads 
to the two conditions, 

5](L,)=0, E("t) = l- 

i ia 

Note that applying this on average in space prohibits us 
from considering charge ordering in our calculations. To 
enforce these constraints we introduce chemical poten- 
tials for the rotors and for the spinous, ^ll and /ij re- 
spectively. This leads to two independent Hamiltonians 



which only talk to each other through the mean fields Xij 
and Bii, 



Ht 



{ij)a 


(15) 


f >; (lUi)-s.. 


(16) 



Ht 






where we've introduced tij which is equal to t on the in- 
plane triangular bonds and equal to t' on the out of plane 
bonds. While the spinon Hamiltonian can be treated 
using mean field theory, the rotor Hamiltonian needs a 
different strategy. 

Two approaches for the rotor Hamiltonian have been 
used in the literature: a self-consistent cluster approactPS 
and a bosonic approachr^. taking the bosonic approach, 
which is most simply tackled using a path-integral for- 
mulation, the imaginary time action takes the form, 

s{e,L) = j dT \yY,Uudre, + ^Li\ 



where fi^ is chosen to eliminate the linear terms in S. 
Due to the symmetry of the action this guarantees that 
(Li) = 0. Next, we integrate out L to get the following 
action, 



s{e) = / dT 

Jo 



2U 



(drt 



^^U.X.,^''^^ 



(18) 
We write this using a bosonic variable 0^ — e*^' subject 
to the constraint that \4>i\^ = 1, giving 

5(0,0; A) - f dr[^^|a,0,|'-*^A, 

^ i i 

+ X! ("^^i^ij ~ 2tijXij) 0i0jj , (19) 

ij 

where A is an auxiliary field introduced to enforce the 
constraint. Treating this new constraint in saddle point 
approximation, the solution of the bosonic part of the 
Hamiltonian is reduced to solving this saddle point equa- 
tion and a free bosonic problem. These saddle point 
equations simply fix the boson number at each site to 
one. 

We assume a uniform state on the Fe sublattice, by 
considering Xij = X ^-nd Bij = B in plane, with Xij = olX 
and Bij = aB out of plane, where a is chosen so that 
we smoothly match the non- interacting limit .I^Sl For the 
Cr spins it is natural to a ssume that the periodicity 
is that found by experiment j^° l ^^l and previous classical 



calculations^, with wavevector of (1/3,1/3,0). Within 
this space of magnetic states, we consider only canted 
classical ground states, that is states where the inplane 
components in each triangle are at 120° but they are 
tilted out of plane by a canting angle ip. These states 
interpolate between a subset of the ground states for 
Jh 3> Jk at ■(/) = and the ferrimagnetic state valid 
for Jk ^ Jh with tp = n/^^. Among these states one 
can see that only those with a finite moment, as one 
sums around a hexagon in the Kagome lattice, will be 
favoured due to the Jk interaction. With these con- 
straints we only have a single set of magnetic states to 
consider parametrized by the canting angle ip. Under 
these assumptions the trimer feels the following effective 
magnetic field (defined as hi — X^aeO "^a) after summing 
around the hexagon 



6Scos^p cos {Q ■ ri)x + ain {Q ■ ri)y 



12S sinipz, 



where Q — (1/3,1/3,0) and S is the magnitude of the 
Cr moment. Note the factor of 2 between the in-plane 
and out-plane components due to the partial cancellation 
as we sum around the hexagon. The phase diagram for 
Jn/t = 0.4 is shown in Fig. [tI and demonstrates the full 
range of phases. 

There are three distinct phases shown in these phase 
diagrams. At low JK/t and U/t we find a metallic phase 
on the Fe sublattice with zero canting angle (i.e. in plane, 
120° ordering) on the Cr sublattice. This metallic state 
is characterized by condensation of the bosonic degree of 
freedom ((0) ^ 0) and a uniform, real Xij — \x\ ™" 
plying the existence of an electron Fermi surface and 
gapless charge excitations. Futhermore, the Fe trimers 
have an induced moment (anti-ferromagnetically) follow- 
ing the Cr moment. This moment scales with JK/t, and 
is thus will be small so long as JK/t is. As we increase 
U/t and keep JK/t small, we find a metal insulator tran- 
sition near Uc/t « 3.5 into a uniform, f/(l) spin liq- 
uid (SL) phase. The metal-insulator boundary is very 
flat, as seen in Fig. [7J , since under the mean fleld 
ansatz we have employed the critical U is only a function 
of X which changes only by a small amount as we in- 
crease Jk (below the jump into the paramagnet phase). 
This SL phase is also characterized by a uniform, real 
Xij but gapped bosons (((/>) = 0), meaning the existence 
of a spinon Fermi surface but gapped charge excitations. 
This insulating phase carries the induced moment as in 
the metallic phase, with the magnitude also being pro- 
portional to JK/t- In both the metallic and insulating 
phases as JK/t is increased (but remains small) the x 
order parameter decreases towards zero. 

From the spin only model of Ref. [20] (accounting for 
the differing normalizations) one expects that for large 
U/t the canting angle will become non-zero at Jk = 2Jh 
and saturate at Jk = ^Jh as one sees in Fig. [7J We find 
that the uniform spin liquid does not support a non-zero 
canting angle. By this we mean that as J/^/i is increased 
the canting angle remains zero until some critical JK/t, 
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FIG. 7: (Color Online) Phase diagram as a function of Jk /t 
and U/t for the value Jh /t ~ 0.4. The inset diagrams show 
a sample of the Cr spin configuration on one of the Kagome 
triangles, x ^^'^ B and the slave-rotor mean field parameters, 
(0) is the rotor condensate and ij) is the Cr canting angle (see 
text). 



wherein the spin liquid is replaced by the canted AF with 
X = B = Q. The converse is not true, as can be seen 
in Fig. [T] We see that the spin liquid phase can be 
destroyed before the onset of the canted magnetic state, 
leaving a window where we have a a simple insulating, in- 
plane antiferromagnet. Once inside this canted AF phase 
with X — B — Q^ the canting angle increases with Jk/^ 
until in becomes saturated and we enter a ferrimagnetic 
phase, with a net magnetic moment. The key feature of 
the phase diagrams we want to emphasize is that for a 
variety of values of Jr and for small Jk this spin liquid 
state is stable and does not coexist with a finite canting. 



V. DISCUSSION 

To discuss the applications of this to FeCrAs the effects 
of fluctuations must be taken into account near the metal 
insulator transtion. These include both charge and gauge 
degrees of freedom and are treated in detail in the work 
of Podolsky et aP^. In this work it is found that there are 
two relevant temperature scales T* and T** that deter- 
mine the qualitative features of the thermodynamic and 
transport properties. These temperature scales vanish 
as one approaches the critical point separating the metal 
and insulator. At temperatures above these scales the 
specific heat has weak logarithmic corrections 



C-Tlnln(l/T), 



(20) 



while the conductivity has a strong temperature depen- 
dence. Specifically, writing cr = ct/ + Ufc where af is the 
spinon conductivity and ui, is boson conductivity, under 
the assumption of weak disorder we have ctj ~ crl™'^ due 



to impurity scattering and 



(21) 



This implies that the effects of the fluctuations on the 
specific heat are much more difficult to discern exper- 
imentally then the effect on the resistivity, which will 
be a monotonically decreasing function of T, as seen in 
the experiments on FeCrAs. This is possible scenario for 
FeCrAs where a nearly linear specific is observed, as these 
logarithmic corrections would only be visible over large 
temperature ranges, where other contributions would be- 
gin to dominate and wash out the signature. This is also 
consistent with the magneto-resistance measurements^, 
where no change is seen in the low temperature resistiv- 
ity under magnetic fields up to 8T, as the magnetic field 
would only couple weakly to the rotor fluctuations. 

A method to test this hypothesis experimentally would 
be a study of the pressure dependence of the trans- 
port and thermodynamic properties. Naively one expects 
that the application of pressure should drive the material 
through the metal-insulator transition, into the quantum 
critical metal and eventually into a Fermi liquid phase. 
Under our scenario, this could in principle be visible as 
the development of a maximum in the resistivity at low 
temperatures as the pressure is increased. Furthermore, 
the specific heat should be relatively unaffected, still only 
being renormalized by a logarithmic term. 

The limitations of the current study deserve some dis- 
cussion as they lead to future directions. Due to the 
natural ansatz used in the slave-rotor study, a number of 
non-trivial spin-liquid states have been excluded from the 
analysis, such as those with a non-trivial phase structure 
or that break translational symmetrjPSl. A full explo- 
ration of the possible phases in this model could yield use- 
ful insights for FeCrAs. Another aspect of this problem 
that requires future work is the inclusion quantum effects 
in the description of Cr, leading to a Kondo-Heisenberg 
model with strong interactions for the conduction elec- 
trons. While the addition of frustrating interactio ns for 
the Kondo spins has attracted attention recently25E3^ 
the inclusion of electronic interactions is largely unad- 
dressed (particularly with frustration on the conduction 
electrons) and is an interesting, but highly non-trivial, 
question for future study. 



VI. SUMMARY 

In this paper we have presented and motivated a min- 
imal model for the low energy degrees of freedom in the 
compound FeCrAs. Starting from the crystal structure 
and using the experimental facts, we have argued that the 
magnetic degrees of freedom are well described by a set 
of classical, localized moments and the electronic degrees 
of freedom take the form of a half-filled Hubbard model 
on the trimer sublattice. The coupling between these two 
subsystems stabilizes a definite magnetic order on the lo- 
calized moments despite the high degree frustration. To 
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(a) xy orbital (b) xz orbital (c) yz orbital 

FIG. 8: (Color Online) The orbitals of the i2 level, rotated 
along the local axes of an As tetrahedron. 



explain the thermodynamic and transport properties of 
this material at low temperatures we propose that the 
electrons residing on the Fe trimers could be close to a 
quantum critical point separating metallic and insulating 
phase. The charge fluctuations associated with the crit- 
ical point strongly renormalize the transport properties 
but provide only small corrections to the thermodynam- 
ics, qualitatively consistent the experimental results on 
FeCrAs. Finally, we have discussed unexplored experi- 
mental consequences of this proposal and future direc- 
tions for theoretical work. 



dron axes. Rotating these into the axes of a tetrahedron 
in a trimer, the orbitals are oriented as shown in Fig. [8] 

From Fig. [8] one can estimate that the only signifi- 
cant hoppings would be between xy — xy, xz — yz and 
yz — xz. One can approach this more quantitatively us- 
ing the ideas of Slater-Koster^S theory to compute the 
orbital overlaps in terms of rotation matrices and irre- 
ducible overlaps. We identify the vector connecting two 
sites in a trimer a as, [x -\- \/3y)/2, giving the (Euler) 
representation R — Rz{0)Ry{—^)Rz{—^). The rotation 
that takes our tetrahedron into the proper axes is given as 
Rrp = jj.(_|)/jg(|)^.(D giving the transformation to 
local axes Rt and Rz{^)Rt for the neighbouring tetra- 
hedron in the trimer. This reduces the number of param- 
eters to three, given by overlaps of ? = 2, m = 0, ±1, ±2 
orbitals displaced along the z direction which we will de- 
note as io-, tjr and ts respectively. The hopping matrix in 
the basis of xz,yz and xy is then given by 



1 


ta-8t„-9t„ 


7t5-16t„-l-9t„ 


^{-7ts-4:t^+3t„ 


32 


7ts~16t^+9t„ 


t^-st^-^gt^ 


^(-7ti+4t„-3t„ 


Vl(7ta+4t„-3t„) 


Vl(-7ti-4t,+3t„) 


i(49i5-12i„+3t<,) 
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Considering the orbitals at atomic separations, we expect 
t(j and ts are positive with t^^ negative. The simplest 
ansatz to try is t^ ^ tg ~ —t^ = t. This gives: 



(Al) 




Appendix A: Hopping integrals 

We first consider direct hopping between the Fe atoms 
in a trimer. The ^2 level is composed of xz, yz and xy 
orbitals with respect to the canonical choice of tctrahe- 



as one might guess by looking at the orbital overlaps in 
the rotated axes. Varying the numerical values for the 
irreducible hopping parameters around this point gives 
qualitatively the same picture as this simple case, justi- 
fying our naive guess. 
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